arXiv:astro-ph/9502075vl 15 Feb 1995 


A HEURISTIC REMARK 
ON THE PERIODIC VARIATION 

IN THE NUMBER OF SOLAR NEUTRINOS 
DETECTED ON EARTH 

H.J. Haubold 

UN Outer Space Office, Vienna International Centre, Vienna, Austria 

and 

A.M. Mathai 

Department of Mathematics and Statistics, McGill University, Montreal, 

Canada 

Abstract. Four operating neutrino observatories confirm the long stand¬ 
ing discrepancy between detected and predicted solar neutrino flux. Among 
these four experiments the Homestake experiment is taking data for almost 
25 years. The reliability of the radiochemical method for detecting solar neu¬ 
trinos has been tested recently by the GALLEX experiment. All efforts to 
solve the solar neutrino problem by improving solar, nuclear, and neutrino 
physics have failed so far. This may also mean that the average solar neutrino 
flux extracted from the four experiments may not be the proper quantity to 
explain the production of neutrinos in the deep interior of the Sun. Occasion¬ 
ally it has been emphasized that the solar neutrino flux may vary over time. 
In this paper we do address relations among specific neutrino fluxes produced 
in the proton-proton chain that are imposed by the coupled systems of non¬ 
linear partial differential equations of solar structure and kinetic equations 
by focusing our attention on a statistical interpretation of selected kinetic 
equation of PPII/PPIII branch reactions of the proton-proton chain. A fresh 
look at the statistical implications for the outcome of kinetic equations for 
nuclear reactions may shed light on recent claims that the ‘ Be -neutrino flux 
of the Sun is suppressed in comparison to the pp- and 8 B neutrino fluxes and 
may hint at that the solar neutrino flux is indeed varying over time as shown 
by the Homestake experiment. 


1 Solar Nuclear Energy Generation: Proton- 
Proton-Chain 

The nuclear energy source in the Sun is believed to be the proton-proton 
chain, in which four protons fuse to form one 4 He nucleus, i.e. 

4p —> 4 He + 2e + + 2v e + Q, (1) 

where Q = M{ 4 He) — 4 M p — 2M e rs 26.73 MeV denotes the energy re¬ 
lease. The three different branches (PPI, PPII, PPIII) to accomplish the 
formation of 4 He in the pp-chain are shown in Figure 1. Neutrinos are 
produced in the pp-chain by nuclear fusion reactions, beta-decay, and elec¬ 
tron capture. The dominant reactions (PPI, 86% of the produced 4 He ) 
produce the great majority of low energy solar neutrinos {$„ SM (pp) « 6.0 x 
10 10 ucm~ 2 s~ 1 ,^ SM (pp) ~ T c -1 ' 2 , where T c denotes the temperature at the 
centre of the Sun) and their number should be a firm prediction of any solar 
model because it is closely tied to the solar luminosity. The second branch of 
the pp-chain (PPII, 14% of the produced 4 He ) yields the ‘Be neutrinos at two 
discrete energies (<3>() SM (jBe ) ~ 4.9 x 10 9 ucm^ 2 s^ 1 , ^ SM (J Be) ~ T c 8 ). In the 
third branch of the pp-chain, the very rare case (PPIII, 0.02% of the produced 
4 He), radioactive 8 B is produced which decays ultimately and is the source 
of high energy 8 B neutrinos (<fA SM ( 8 B) tv 5.5 x 10 6 z/cm -2 s -1 , ^^ SM ( 8 B) ~ 
T 48 ). The specific neutrino fluxes are taken from the Standard Solar Model 
(SSM) of Bahcall and Pinsonneault (1992). The overall energy produc¬ 
tion of the pp-chain is Q = 26.73 MeV, however, the three branches pro¬ 
duce a different amount of energy (PPI: E 7 = (26.73 — 0.265)MeV, PPII: 
E-y = (26.73 — 0.861)MeV, PPIII: = (26.73 — 7 )MeV) due to the energy 

loss carried off by elusive neutrinos. The fluxes of high energy solar neutrinos 
are especially sensitive to the central temperature, T c , mainly because of the 
energy dependence of the cross sections of the respective nuclear reactions. 
The branching ratios (the percentage of which each branch of the pp-chain 
contributes to the production of 4 He) are strongly dependent on the nuclear 
reaction probabilities and on the density and temperature profiles inside the 
Sun. Assuming that the Sun is in a state of quasistatic equilibrium, the solar 
luminosity L Q tells us the total energy generation rate which can be turned 
into a constraint on the total solar neutrino fluxes, that is 

Lq = 13.1 ((j>u{pp) - 4>v( 7 Be) - 0v( 8 7?)) + 25.6<%( 7 .Be) + 19.5<%( 8 .B). (2) 



The luminosity L 0 observed at the current stage of evolution of the Sun 
corresponds to the energy that was generated in the gravitationally stabilized 
solar fusion reactor 10 7 yr ago (Helmholtz-Kelvin timescale). The quasistatic 
assumption allows one to equate the present luminosity with the present 
nuclear energy production rate (Mathai and Haubold, 1988). Normalizing the 
neutrino fluxes in (2) to those of the Standard Solar Model ($„ = <Pv/(pv SM ) 
leads to a luminosity constraint indicating the degree of contribution of the 
respective neutrino flux to the total solar neutrino emission: 

1 = 0.913$„(pp) + 0.071& u ( 7 Be) + 0.00004^( 8 5). (3) 

To reveal how the pp-chain operates in the Sun is to measure the individ¬ 
ual neutrino fluxes <pu(pp), (pv('Be), and (p u ( 8 B) in (3), thereby fixing the 
branching ratios of PPI, PPI1, and PPI11 as indicated in Figure 1. Four 
experiments are in operation now to accomplish this solar neutrino spec¬ 
troscopy. Kamiokande measures exclusively the high-energy flux (pvi^B ), 
thus PPI11, in real-time spectroscopy. Homestake observes primarily <p u ( 8 B) 
and to a much lesser extend (p v { 7 Be ), that is, the branching of PPI1 and 
PPI11. GransSasso/Baksan detects primarily the low-energy flux (p u {pp) and 
to a lesser extend the fluxes (p(‘ Be) and (p^( 8 B ), thus focusing on the branch¬ 
ing of PPI and PPII/PPIII. 

2 Spatial Distribution of Solar Neutrino Sources: 
Standard Solar Model 

Figure 2 shows the neutrino production as a function of the dimensionless 
distance variable x = R/R Q , starting from the center of the Sun for the 
Standard Solar Model (Bahcall and Pinsonneault, 1992). We note that ex¬ 
periments for the detection of solar neutrinos are looking into different depths 
of the solar core as they are only sensitive to specific neutrinos produced in 
the respective nuclear reactions. The region in which the low-energy pp- 
neutrino flux is produced is very similar to that of the total nuclear energy 
generation. Because of its strong temperature dependence, the high-energy 
8 B-neutrino production is peaked at the very small radius x = 0.05 and is 
generated in a much narrower region in comparison to the other two neutrino 
sources (Table 3). 



The contribution of the neutrino fluxes <f> u (pp), <f>v{ 7 Be), and <f) u ( 8 B) to 
the Homestake, Kamiokande, and GranSasso/Baksan experiments is 

^(Homestake) ~ 6.20 I/ ( 8 i?) + 1.2 Be)SNU, (4) 

Kamiokande) = 4> U ( 8 B)^ I/ ( 8 B), (5) 

^(GranSasso/Baksan) ~ 13.80^( 8 il) + 35.8 (f) u (‘ Be) ± 70.8 4> u (pp)- (6) 

The coefficients in Equation (4)-(6) are corresponding to the capture rates 
predicted by the Standard Solar Model for each respective solar neutrino 
source. 

Table 1 summarizes the predictions of the total capture rates of the Stan¬ 
dard Solar Model for the Kamiokande, Homestake, and GranSasso/Baksan 
experiments. The Kamiokande solar neutrino flux is given in units of 
10 6 wm“ 2 s _1 , while the Homestake and GranSasso/Baksan rates are given in 
solar neutrino units (1SNU = 10 _36 ^afom _1 s _1 ). The uncertainties shown 
in Table 1 are lcr. 


Experiment 

Kamiokande 

Homestake 

Gr anSasso/Baksan 

Predicted 
capture rate 

5.69 ±0.82 

8.0 ±3.0 

131.5±ft 


Table 1. 


3 Four Solar Neutrino Experiments: Solar 
Neutrino Problem and Be-Neutrino Defi¬ 
ciency 


Experiment 

Reaction 

Energy 

threshold 

Location 

Kamiokande 

v e ± e~ —> v e ± e~ 

7.5 MeV 

Japan 

Homestake 

v e ± 37 Cl —U 7 Ar ± e~ 

0.814MeU 

USA 

GranSasso 

(GALLEX) 

v e ± 71 Ga —> 71 Ga ± e~ 

0.233 MeV 

Italy 

Baksan 

(SAGE) 

v e ± 71 Ga —A 1 Ge ± e~ 

0.233 MeV 

Russia 
























Table 2. 

There are four experiments currently operating to detect neutrinos com¬ 
ing from the Sun (Table 2). The Kamiokande experiment is a water Cerenkov 
detector which measures the energy of the scattered electrons (Nakamura, 
1993). Due to its energy threshold it is only sensitive to the high-energy 8 fi 
neutrinos from branch PP111 of the pp-chain. The Homestake experiment 
consists of 10 5 gallons of C- 2 CU and detects solar neutrinos via capture on 
the chlorine (Davis, 1993). Its energy threshold allows to detect the higher 
energy line of 'fie-neutrinos from branch PPII as well as the high-energy 8 fi 
neutrinos from branch PPIII. The two gallium experiments at GranSasso and 
Baksan are sensitive to the low energy pp-neutrinos from the PPI branch of 
the pp-chain as well as to the higher energy 'fie- and 8 fi-neutrinos (Ansel- 
mann et ah, 1994; Abdurashitov et ah, 1994). The predicted contributions to 
the Homestake and GranSasso/Baksan experiments based on the Standard 
Solar Model are shown in Table 3. 


Neutrino 

source 

Homestake 

experiment 

Percentage 
of total 
capture rate 

GranSasso/Baksan 
experiments 

PP 

0.0 

= 0% 

70.8 

pep 

0.2 

- 2.5% 

3.1 

7 fie 

1.2 

= 15% 

35.8 

8 fi 

6.2 

77.5% 

13.8 

i8 N 

0.1 

:z 1.25% 

3.0 

15 0 

0.3 

3.75% 

4.9 

total 

capture rate 

8.0 ±3.0 

= 100% 

131.5±i7 


Table 3. 

The experimental results of the four solar neutrino experiments are given 
in Table 4 and can be compared with the predictions of the Standard Solar 
Model as shown in Table 1. 




Experiment 

Detected 
capture rate 

Detected capture rate/ 
predicted capture rate 

Kamiokande 

2.89^ ±0.35 

0.50 ±0.07 

Homestake 

2.55 ±0.17 ±0.18 

0.32 ±0.03 

GranSasso/ 

Baksan 

77 ±9 

0.59 ±0.07 


Table 4. 

From Table 4 it is evident that the results of the four experiments are between 
1/3 and 1/2 of the neutrino capture rates predicted by the Standard Solar 
Model. This deficit of solar neutrinos is called the solar neutrino problem 
which poses a serious conflict with the constraint of the overall solar luminos¬ 
ity in (2) and (3). Additionally, the comparison of the three detected capture 
rates in Table 4 with the predicted capture rates in Tables 1 and 3 shows that 
the Kamiokande rate is less suppressed than the Homstake rate. Because the 
Homestake experiment has a lower energy threshold, the lower detected cap¬ 
ture rate suggests that the 7 B e-neutrinos are more suppressed than the high 
energy 8 Z?-neutrinos. However, any reduction of the ‘Be production rate by 
lowering the temperature T c would affect immediately both the 7 Be and 8 B 
neutrino production equally. This fact seems to pose an additional problem 
in finding a solution of the solar neutrino problem in terms of solar, nuclear, 
and neutrino physics on which the Standard Solar Model is based. This is 
particularly true for the so-called cooler Sun models (Castellani et al., 1994). 


4 Argon-Production Rate of the Homestake 
Experiment: Variations Over Time 

Figure 3 shows the 37 Ar production rate detected by the Homestake exper¬ 
iment from 1970.8 to 1991.6 (Davis, 1993). The average 37 Ar production 
rate (combined likelihood function) for the 94 individual runs shown was 
0.509 ± 0.031 argon atoms per day. Subtracting a total background 37 Ar 
production rate of 0.08 ± 0.03 atoms per day yields the production rate that 
can be ascribed to solar neutrinos: 0.429 ±0.043 atoms per day or 2.28 ±0.23 
SNU (the rate in SNU is equal to 5.31 times the captures per day in the Home- 
stake experiment). This average capture rate is commonly compared to the 




predictions of the Standard Solar Model as shown in Tables 1 and 3 for the 
pp-chain and the CNO cycle. This procedure does not take into account in 
any way the apparent time variation in the observed 37 Ar production rate 
evident in Figure 3. 

Figure 4 shows a five-point moving average of the 37 Ar production rate, 
removing high frequency noise from the actual time series collected in the 
Homestake experiment as shown in Figure 3. One notes in the five-point 
moving average that in the periods 1978 to 1979 and 1987 to 1988 a su- 
pression of the 37 Ar production rate seems to occur. The overall shape of 
the five-point moving average suggests that there are two distinctive epochs 
spanning the time periods 1971 to 1980 and 1980 to 1989. Each epoch shows 
a shock-like rise and subsequent rapid decline of the 37 Ar production rate. 
Further, the five-point moving average of the 37 Ar production rate reveals 
that each of the two distinct cycles covers a time period of around 9 years. 
Each cycle exhibits a slow but shock-like increase, reaching a peak, succeeded 
by a rapid decrease to a minimum value of the 37 Ar production rate. This 
pattern is repeated for a second nine-year period and seems to start for a 
third period in 1989 (Haubold and Mathai, 1994). Each of these cycles can 
be reproduced by a mechanism discussed in the following Section. Fourier 
analysis of the 37 Ar production rate data in Figure 3 reveals a power spec¬ 
trum showing the harmonic content in this time series in terms of a series 
of distinctive periodicities which is shown in Figure 5. Fourier analysis also 
indicates that the harmonic content in the 37 Ar production rate data is dom¬ 
inated by periodicities of 0.57, 2.2, 4.8, and 8.3 years (Haubold and Gerth, 
1990). 

5 Kinetic Equations: Lifetime Densities 

The production and destruction of nuclei in the proton-proton chain of re¬ 
actions can be described by kinetic equations governing the change of the 
number density N t of species i over time, that is, 

iNi = - N i N i < av > ij + N k N l < OV >kl, ( 7 ) 

dt j k,l^i 

where < av > mn denotes the reaction probability for an interaction involving 
species m and n, and the summation is taken over all reactions which either 



produce or destroy the species i. The first sum in (7) can also be written as 
- Y N i N j < av >ij= ~ N i(Y.] N J < av >ij ) = ~ N i a U ( 8 ) 

3 3 

where a* is the statistical expected number of reactions per unit volume 
per unit time destroying the species i. The reciprocal of the quantity a* is 
the lifetime of species i for interaction with species j for all j. It is also a 
measure of the speed in which the reaction proceeds. If the reaction results 
in the production of a neutrino, for example, then the reciprocal of cq is the 
expected time it takes to produce this neutrino in the solar interior. In the 
following we are assuming that there are Nj(j = 1,of species j per 
unit volume and that for a fixed N t the numbers of other reacting species 
that react with the i-th species are constants in a unit volume. Following the 
same argument we have for the second sum in (7) accordingly, 


+ Y N k N i < av >ki= +NA, 

k,l^i 

where Nfi, is the statistical expected number of the i-th species produced per 
unit volume per unit time for a fixed N t . Note that by nature the number 
density of species i, N t = N t (t), is a function of time while the < av > mn 
are supposed to depend only on the temperature but not on the time t and 
number densities N r Then equation (7) implies that 

j t N i {t) = -{a i -b i )N i {t). (9) 

For equation (9) we have three cases, q = a* — 6* > 0, c t < 0, c x — 0, of 
which the last case says that N t (t) does not vary over time, which means 
that the forward and reverse reactions involving species % are in equilibrium. 
The first two cases exhibit that either the destruction (c t > 0) of species i or 
production (c t < 0) of species i dominates. 

For the case q > 0 we have 


d_ 

dt 


Ni{t) 


-CiNiit ), 


and it follows that 


Ni(t)dt = Ni(0)e~ Cit dt, 


( 10 ) 



where iVj(O) is the number density of species i at time t — 0. If c, in (10) is 
a function of time, say c*(t), then c,t in (10) is to be replaced by / dtci(t). If 
the arrival distributions for the other species are Poisson, then q(f) will be 
of the form df, where dj > 0 independent of t. In this case the exponent in 
(10) is / dttCiit ) = dit 2 / 2. Contrarily, when c % is a constant, the total number 
of reactions in the time interval 0 < t < f 0 is given by 


/■to rto NAO) j. 

/ dtNi(t) = Ni(0) dte~ Cit = lK \l e~ Cit0 ). 

Jo Jo Ci 

(ii) 

In (11), 1 — e~ Cit ° is the probability that the lifetime of species i 
t has the density 

is < t 0 when 

fit) = Cie~ Cit , 0 < t < oo, ^ > 0, 

(12) 

or 

n&) = 

Ci 

(13) 

When Ci = q(I) = d,t then 


m) = [A ) 1 ' 2 

(14) 

where ^ 

hit) = f — -'j e ~dit 2 / 2 ^ o < t < oo, di > 0. 

(15) 


The density in (12) will be called the lifetime density for the destruction of 
species i, with the expected mean lifetime 


E(t) = 1. (16) 

If the lifetime density is as given in (15) then 

/ 2 \'/ 2 

B w - U ■ < 17 > 

From (12) and (16) we can make the following observations: 

(i) q can be interpreted as a measure of net destruction, the larger 

the value of c t the faster the net destruction. 



(ii) Ni ^ f(t) At can be interpreted as the amount of net destruction over 
the small interval of time At. The faster the net destruction the shorter the 
lifetime. 

(iii) The quantity 

r dt N ‘(°) _ w) 

Jo 


-m = 


(18) 


Q Ci 

can be interpreted as the total net destruction of species i starting with the 
initial number iVj(O). 

(iv) If the net destruction of species i produces a species k, for example 
a neutrino, then the number produced is proportional to Ni ^ . 

If the lifetime for the production of a species k due to the net destruction 
of species i is denoted by r, then r is a constant multiple of t, say r = ait, 
where t has the lifetime density f(t). But the densities of t and an t (an > 0) 
belong to the same family of distributions and hence the density of r can be 
written as 

f( r ) — Qi e ~ 6iT i 0i > 0,r > 0, (19) 

where 0, = Cj/ai and thus the total production is 


6 Dampening of Reactions: Poisson Arrivals 

Suppose that after a certain period of time of net destruction, say t 0 , a 
dampening effect starts to slow down the net destruction of species i with 
initial number iVj(O). Let this dampening variable be denoted by t 2 , where 
T 2 is again proportional to the lifetime, say a 2 f. Then the lifetime density 
associated with r 2 is of the exponential type, belonging to the same family as 
in (19). Let T\ and r 2 be independently acting or statistically independent. 
Let the delay in time for r 2 to start be c = a 2 f 0 and let the densities of T\ 
and r 2 be denoted by 

fArj) = fa ljTj ■ ~i > 0,/3j > 0 ,j = 1,2 (20) 

where [A — 9 t — Q/ai of (19) and let /3 2 = Cj/a 2 . Then the net destruction 
of species i is proportional to u = T\ — (r 2 — c) = T\ — r 2 + c with the joint 
desity of T\ and r 2 given by 


f(n ,t 2 ) = (3i(3 2 e~^ 1T1+ ^ 


1 Jj > 0 ,(3j > 0 ,j = 1,2 


( 21 ) 



due to the statistical independence of T\ and t 2 . The density of u, denoted 
by g(u), is the following (Mathai, 1993) 


9 0 ) 


0102 
01+02 


e c \ c < u < 00 


0102 

01+02 


e 02(u c) j _ QO < u < c 


( 22 ) 


where 


P1P2 = Cj 
/3\ + /3 2 Oi + o; 2 


observing that /?i = q/cki and /3 2 = Cj/a 2 . If the net destruction of species 
i is exceeding the dampening rate, then > (3 2 and the following Figure 6 
illustrates the behaviour of the density u in (22). 

Figure 6 shows a non-symmetric Laplacian, slowly rising and rapidly falling. 
At the time t = t 0 the net destruction of species i is given by Ni ^ 0 ) (1 — e~ Cito ). 
Then the production of species, for example neutrinos, as a result of the net 
destruction of species i, in an instant of time is given by 


aiQ ! 2 



(23) 


which is a constant multiple of g(u), where g(u) is given in (22). Hence 
the shape of the curve for the net destruction of species i and the resulting 
production of species k will be the same as of g(u) shown in Figure 6. The 
production of resulting species in a small intervall of time At is Ag(u)At with 
to = c/ck 2 starting with a constant initial number iVj(O) of species i, where 


A — (X\(X2 


Nj( 0 ) 

Ci 



(24) 


since dug(u ) = 1. Here the integration is done from — 00 to + 00 . Note 

however, that when c is large enough the probability for u being negative 
will be negligibly small and hence A in (24) is a good approximation to the 
total production. 

We observe in (24) that when c is small, A is small and A is an increasing 
function of c as shown in Figure 7. 

Note that A in (24) is the result of assuming that the initial number 
iVj(O) of species i per unit volume is a constant. If the species i is arriving 



to the unit volume according to a Poisson distribution with parameter A* 
(Poisson arrivals), then iVj(O) in (24) as well as in the previous formulae is 
to be replaced by its expected value, that is i?[iVj(0)] = \ in the considered 
case. In Poisson arrivals one can take the expected number to be \ = 7 {t, 
where t is the duration of destruction and 7 ; is a constant independent of 
time t. In this case A in (24) becomes 


A 


Qi\OL 2 


CJ 


77(1 


(25) 


where to — c/a 2 is the time where the dampening effect starts. 


7 Proton-Proton Chain: Branches II and III 


The fusion of four protons to produce one helium nucleus in the pp-chain 
is accomplished in at least three different branches in the chain (Figure 1). 
This branching results in uncertainties of the predictions of the 'Be- and 8 B 
neutrino fluxes in the Standard Solar Model and needs particular attention in 
discussing the results of those solar neutrino experiments which are detecting 
exclusively ' Be- and S B neutrinos (Homestake and Kamiokande experiments 
in Tables 1,2, and 4). Without any branching in the pp-chain, the number 
of all reactions and neutrinos would be equal, that means 

4>v{pp) = = N. 


As shown in Figure 1, 3 He can interact with another 3 He nucleus to produce 
right away 4 He (PPI branch), or 3 He can fuse with 4 He to produce a 'Be 
nucleus and subsequently to open branches II and III of the pp-chain. The 
branching ratio r is determined by the reaction probabilities < av >ij and 
number densities N t : 


< CTV >34 


A4 

'AV 


(26) 


1 - r < av >33 

where the notations have been explained in the preceding section. With 
regard to branches II and III in Figure 1, 'Be can capture an electron to emit a 
'Be neutrino, or it can fuse with a proton to produce 8 B which immediately 
decays and produces a 8 B neutrino. The branching ratio r' for PPII and 


PPIII is 

r' < av >17 Ni 
< av > e7 ‘ N e ' 


1 — r' 


( 27 ) 



With (26) and (27) the following relations between the number of chains 
to produce 4 He and the neutrino fluxes produced by the three branches are 
established, 

N 

Mpp) = y (2 - r )> ( PPI ), (28) 

<j )u CBe) = jr(2-r'), (PPII ), (29) 

M 8 B) = ^rr', (PPIII ). (30) 

Equations (28) to (30) show the link of the three branches of the pp-chain 
which is eventually governed by the reaction probabilities < av > ^ and num¬ 
ber densities N t in the system of kinetic equations in (7) and by the profiles 
of density and temperature of the solar model. Basic assumptions for equa¬ 
tions (28) to (30) are that the Sun is in thermal equilibrium which fixes the 
number of chains through (2) and that the nuclei responsible for neutrino 
production are in thermal equilibrium with the ambient plasma which al¬ 
lows to determine the neutrino fluxes by the reaction probabilities. For the 
latter assumption the characteristic time for significant energy exchange by 
Coulomb collisions between reacting species must be orders of magnitude 
less than the characteristic time it takes to produce a neutrino in the solar 
interior (Maxwell-Boltzmann reaction rates). These basic assumptions still 
leave the question open on what is relevant for branching governed by kinetic 
equations: The time for reducing the protons to thermal equilibrium with 
the ambient plasma (~ 10 -20 yr) or the lifetime of a proton to undergo a 
reaction with a second proton to produce, among other species, a neutrino 
(cs 10 10 yr)? This question will be addressed in the following section for three 
reactions of the branches II and III of the pp-chain. 


8 Production - Dampening Mechanism: Lapla- 
cian Behaviour 

Consider three sets of Laplacians of the type given in Figure 6, one set 
consisting of one Laplacian with to = |(1) units of time, five successive 
Laplacians with to = |(0.2) units of time each in the second set, and the 
third set consisting of eight successive Laplacians with to = \ (0.125) units 



of time each. Suppose that we consider the Laplacians for a total arbitrary 
time interval of t — 1 unit of time. Let the total destruction of species i by 
one Laplacian of set 1, the five Laplacians of set 2, and the 8 Laplacians of 
set 3 be denoted by Ai,A 2 , A 3 respectively. Then we have from (25) 


A 1 = 


OL\OL2 2 


7?[1(1 ) 2 ] (i - e-^ (1) ) , 


•V = ^7 2 [5(0.2)=] (1 


A* — 


c 2 h 
CX\OL 2 2 


g-Ci^O.2) 


(31) 


^7 2 [8(0.125) 2 ] (l - . 


If c* is large so that e Ci ^ is negligible, then the total contributions coming 
from the three sets are respectively, 


Ai 


Ai + A 2 + A% 


= 0.755, 0.15, 0.095, j = 1,2,3 


(32) 


respectively, that is, 75.5%, 15%, and 9.5% for each of the three reactions. 

The Laplacians can be approximated by using triangles in the following 
way. From Figure 6 it is noted that the maximum height of the Laplacian is 
at u = c which will then be 


Pi P2 _ Ci 
Pi + P2 ex 1 + a 2 


(33) 


Suppose that a 2 = 3cki, which will imply that Pi — ^ and P 2 = 3 ^-, which 
in turn means that the net destruction rate is three times the dampening 
rate. Suppose that Pi = \/3 b, where b > 0 is a constant and t 0 = | b. Then 
the maximum hight of the Laplacian in Figure 6 is 


P1P2 = Cj = V 3 b 
Pi + P2 CKi + Q?2 4 

In this case the Laplacian approximates to the following triangle shown in 
Figure 8 . 

If we take three sets of triangles shown in Figure 8 , where the first set consists 
only of one triangle with 6=1 time unit, the second set contains 5 successive 
triangles with 6 = 0.2 time units each, and in the third set there are 8 
successive triangles with 6 = 0.125 units each, and if the total areas of these 



three sets of triangles are denoted by Ai, A 2 , and A 3 similiar to (31), then the 
respective areas are in the proportion 75.5, 15, and 9.5 percent respectively 
(see Table 5). 


Total area 

Contribution to the 

9 year cycle 

3/ Ar production rate of 
the Homestake experiment 
(Table 3) 

1 triangle 
Reaction 1 

75.5 % 

8 B contributes 77.5% 

5 triangles 
Reaction 2 

15.0% 

'Be contributes 15% 

8 triangles 
Reaction 3 

9.5% 



Table 5. 

Triangles as the ones shown in Figure 8 have been used in the following 
graph, where in Figure 9, a and f3 denote the shifts in the starting point of 
the set of 5 and of the set of 8 triangles, respectively. The starting point of 
the one big triangle has been chosen as t — 0. 


9 Conclusion 

The time variation of the argon production in the Homestake experiment 
which is ascribed to be produced by solar neutrinos can be explained as 
follows. The original Homestake data, but more distinctively the five-point 
moving average of the data, seem to show cycles of approximately nine years 
duration. The reactions of the PPI1 and PP111 branches of the proton- 
proton chain are producing neutrinos through the 7 Be- and 8 B- reactions. 
If one assumes that a dampening mechanism operates for three reactions 
of the PP11 and PP111 branches as discussed above, the variations of the 
argon production in the Homestake experiment over time can be explained 
on purely statistical arguments based on lifetimes and their ratios for the 
three reactions. For these nuclear reactions, destruction and dampening may 
work opposite to each other. If the 

destruction rate is approximately three times the dampening rate, if the 
destruction is \/3b for some b > 0, and if the dampening effect starts |b 




time units from the starting time t — 0, then the time variation cycles seen 
in the argon production in the Homestake experiment can be reproduced 
by considering a scenario of three sets of reactions of the PPII and PPIII 
branches of the proton-proton chain, one set with 6 = 1 unit of time, say 
9 years, the second set consisting of 5 successive reactions with b = 0.2 
time units each, and the third set consisting of 8 successive reactions with 
b = 0.125 time units each. 
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Table captions: 


Table 1 

Table 2 
Table 3 

Table 

Table 5 


: Predictions of the Standard Solar Model for the Kamiokande, 
Homestake, and GranSasso/Baksan experiments (Bahcall and 
Pinsonneault, 1992). 

: The four currently operating solar neutrino experiments. 

: Predicted capture rates in SNU from various flux components of the 
pp-chain and CNO cycle for the Homstake and GranSasso/Baksan 
experiments. The uncertainties are the total theoretical range, 

~ 3<r (Bahcall and Pinsonneault, 1992). 

: Comparison of the detected rates of the four solar neutrino experiments 
with the predicted rates of the Standard Solar Model (Bahcall and 
Pinsonneault, 1992). The Kamiokande flux is in units of 10 6 z/cm _2 s _1 , 
while the Homestake and GranSasso/Baksan rates are in SNU. 

: For the destruction-dampening mechanism considered here, the area 
of the triangle governing reaction 1, the combined areas of the 
5 triangles for reaction 2, and the combined areas of the 
8 triangles of reaction 3, are proportional to the combination 
of three sources to the total capture rate of the Homestake experiment. 



Figure captions: 

Fig. 1: The proton-proton chain and its three different branches to accomplish 
to formation of A He. 

Fig. 2: The neutrino production as a function of the dimensionless distance 
variable x = R/Rq in the Standard Solar Model of Bahcall and 
Pinsonneault (1992). 

Fig. 3: The argon-production rate detected by the Homestake experiment from 
1970.8 to 1991.6 (Davis, 1993). 

Fig. 4: The five-point moving average of the argon-production rate data as 
shown in Figure 3. 

Fig. 5: Power spectrum of the argon-production rate data in Figure 3 obtained 
by Fourier analysis. 

Fig. 6: A non-symmetric Laplacian, slowly rising and rapidly falling. 

The function g(u) is the density of u, describing the destruction- 
dampening mechanism for nuclear reactions involving species i. 

Fig. 7: The behaviour of the function A in (24) denoting the total destruction of 
species i. 

Fig. 8: The non-symmetric Laplacian shown in Figure 6 can be approximated by 
a non-symmetrical triangle. 

Fig. 9: The destruction-dampening mechanism in (25) and (26), where the Laplacians 
has been approximated by triangles, can reproduce the time variation of 
the original argon-production rate within the range of the error bars 
attached to each run. 
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